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Abstract 

We present a numerical determination of the order-disorder interface tension, a Q d, for 
the two-dimensional seven-state Potts model. We find a d = 0.0114 ± 0.0012, in good 
agreement with expectations based on the conjecture of perfect wetting. We take into account 
systematic effects on the technique of our choice: the histogram method. Our measurements 
are performed on rectangular lattices, so that the histograms contain identifiable plateaus. 
The lattice sizes are chosen to be large compared to the physical correlation length. Capillary 
wave corrections are applied to our measurements on finite systems. 



Phase transitions of interest to cosmology and colliders may well turn out to be of first 
order. For this reason there has been increased interest in first-order transitions. Among 
the parameters governing the dynamics of such transitions are the interface tensions. These 
are properties of the equilibrium system. There have been attempts to measure these on 
the gauge theories of relevance to high-energy physics. Recently, however, the control of 
systematic biases in these measurements have been questioned. Such uncertainties call for 
detailed tests of the methodology in models which are more tractable than gauge theories. 

The two-dimensional Potts models are used very often for such purposes. They are also 
interesting because of connections with exactly solved two-dimensional models in statistical 
mechanics. They are defined by the partition function 

Z = Trexp(-/3F), where H = -£<W 3 - (1) 

The spins Oi at each site i of a lattice can take q different values, and the notation (ij) 
denotes nearest neighbour pairs. There is an order- disorder phase transition in the models 
for all q, with a singlet magnetisation as an order-parameter. This transition occurs at a 
coupling f3 c = \og(y/q + 1). For q > 4 the transition is of first-order. 

The coexisting phases can be characterised by the internal energy densities 

E = y(H), (2) 

where the angular brackets denote thermal expectation values and V is the volume of the 
lattice. In the ordered and disordered phases at /3 C , the values of the internal energy densities 
will be denoted by E° and E d respectively. For the Potts models these values are known 
exactly [3J. Recall that as (3 — ► we obtain a pure ordered phase with E = E°. As the 
system is heated, its internal energy increases at constant f3. This is achieved by creating 
larger and larger volume fractions of the disordered phase. Coexistence of the two phases 
implies the creation of interfaces. Finally, when E = E d , a purely disordered phase is 
obtained and with additional heating E increases with a decrease of (3. In the coexistence 
region, E° < E < E d at j3 c , there is a non-extensive part of the free energy. This part can 
be used to define the interface tension. 

A recent exact computation of the correlation length at f3 c in all Potts models with q > 4 
has been identified with the correlation length in the disordered phase, on the basis 
of a large-g expansion ||. A duality argument || then gives the interface tension between 
two different ordered phases, a QO , as 

Ooo = I /id- (3) 
Now making the assumption of perfect wetting, 

a OQ = 2a od , (4) 



one obtains the prediction 

a od = 1/2&. (5) 

Since ^ is exactly known, a measurement of a Q d would check the relation in eq. ([5]) and thus 
test the conjecture of perfect wetting in eq. (|J). This assumption is for large q by the proof 
of the opposite inequalities <j 00 > 2a Q d and <j 00 < 2a od |J. A recent measurement at q = 15 
H has used the argument outlined here to verify that perfect wetting holds for that model. 

The status of the model with q = 7 is interesting. With the assumption of perfect wetting, 
the exactly known value of ^ would predict a od = 0.010396 • • •. A first measurement using 
the so-called integral method gave o Q d = 0.0943(6) J7j. A value in agreement, a od ~ 0.10 
was obtained by a different technique. The surprisingly large violation of perfect wetting 
led to criticism of these two techniques for the measurement of interface tensions. Two 
measurements using the histogram technique were then performed. The first gave a Q d = 
0.0121(5) fH and the second a od = 0.01174(19) flOfl . Although they are in agreement with 
each other, they are in violation of perfect wetting by about seven standard deviations. The 
lattice sizes used in both these studies were not much larger than and hence the results 



may require correction. This was also noticed in [ XT | . In this letter we report a study of the 
q = 7 model by the histogram technique correcting for all sources of systematic errors. 

Consider a finite lattice of size L x L z (L z > L), with periodic boundary conditions 
imposed on the configuations. The histogram of the probability distribution of E, Pl(E)cIE, 
was first used for the extraction of an interface tension in ||12j| . Pl{E) has peaks at E° and 
E d , shifted through finite size effects. In an interval of E bounded by these two values, there 
is a minimum, p™ n ; and a nearly flat region around it. This region is due to two interfaces 
spanning the smaller direction, L, and seperating an ordered from a disordered phase. The 
histogram method consists of a model for p™ n and the nearly flat plateau contained in 
the interval E° < E < E d . Such a model contains as a parameter the interface tension of 
interest, which is then extracted by comparison with data. 

In the region of interest, Pl{E) is independent of E as long as the free energy is inde- 
pendent of the relative positions of the two interfaces. Fluctuations of the interface shapes 
are described by the capillary wave model . For non-interacting interfaces, one obtains 



H 



Pf n = ^exp(-2a od L). (6) 

This contribution takes into account interfaces perpendicular to the z-direction. For L z > L, 
interfaces parallel to the z-direction give small corrections to this result. For the lattices we 



have chosen to work with, these corrections are negligible. In [[HJ it was shown that the 
formula in eq. (|5|) holds for two-dimensional systems even when the global condition that E 
lies between the two pure-phase energies is imposed. 



Whenever all the above conditions are fulfilled, one can extract the interface tension from 
either of the two quantities defined below. The pre-exponential factors in eq. cancel in 



1 , 31nL 31nL 2 



F} = —\n-!—--— + --^^. (7) 

where pj™ 1 ^ [ s the geometric mean of the maxima of Pl(E) corresponding to the ordered 
and disordered phases. Taking into account capillary wave fluctuations of the interfaces, we 
obtain the finite size scaling formula 

Fl = Vod + j for L z >L^oo. (8) 

As a cross-check one can also measure a second quantity from the histograms, 

Capillary wave fluctuations of the two interfaces also yield the finite size scaling formula for 
this quantity — 

Fl = °od + j for L z > L -> oo. (10) 

We shall use both the quantities, F[ and Fl, along with their finite size scaling formulae to 
estimate the interface tension. 
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1.2931 
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140 x 280 


1.293432 


20.0 




1.293562 


12.0 



Table 1: The lattice sizes, couplings and the number of sweeps, N, used in this study. 



Care is required to ensure that the relations above may be applied to the data obtained. 



A sufficient check is to confirm that Pl(E) indeed has a flat plateau at its minimum. In [16 
it was found that the condition L z > L was necessary to obtain non-interacting interfaces. 
Previous studies of this model had used only square lattices. In this study we consistently 
use the value L z = 2L. It has turned out that the lattices need not be elongated further for 
a flat plateau in the histogram to develop. A second necessary condition is that the system 
break up into several domains. This requires L z ^> £. All our lattices satisfy this condition, 
whereas earlier studies were performed on much smaller lattices. Finally, we note that the 
capillary wave model of interface fluctuations has been applied successfully in other contexts, 
and is bound to be an improvement over the neglect of interface fluctuations. 




Figure 1: Measured probability distributions for the energy density E. The dotted line shows 



the result of a computation following [|TTJ for the largest lattice 



The simulations were performed with L x 2L lattices. Since the longest physical corre- 
lation length at (3 C is = 48.09 • • •, the minimum value of L used was 60. For L < 100 
we used a cluster update. For larger values of L the performance of this algorithm was not 
satisfactory. Hence we used a multicanonical version of a two hit Metropolis algorithm in 



the form specified in [[L7[]. The energy was measured every tenth sweep. Between 1000 and 
5000 initial configurations were discarded for thermalisation. The couplings were chosen so 
that the canonical weights were roughly equal for the two phases. The second run on the 
L = 140 lattice was an exception, since this was performed at (3 C . The data was reweighted 
to the canonical distribution at f3 c for each lattice. A cross-check on the statistics of our runs 
was that the weights of the two peaks at f3 c are in the ratio 1:7. For the L = 140 lattice, the 
results presented later average over the two runs. 

We show the histograms obtained in our computations in Figure |I| We would like to 
emphasise that flat plateaus are obtained over an interval between the two peaks for L > 80. 
In a criticism of the histogram method as applied to the ten-state Potts model on square 



lattices |pT| , it was argued that Pl(E) could be reproduced by an analytic computation 
which neglects all interfaces which span the lattice. This is not true for the lattice sizes and 
shapes we choose to work with. In Figure [I] we show our histograms and a comparison with 



a computation following the methods of [ITT | . The analytic computation strongly underes- 
timates the probability in the intermediate region of E. We also find that the discrepancy 
increases with increasing L. This is evidence that interface contributions are important in 
this region of E. 
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Figure 2: Results for F l L together with linear fits. The crosses are results for square lattices 



from P, pi 



We extract Fl and F% from the histograms. A jack-knife procedure is used for the 
estimates of the averages and errors. The results quoted were obtained using between 6 and 
24 jack-knife blocks. We checked that the results were stable under change in the number of 
blocks. A finite size scaling analysis of each of these measurements is performed by fits to 
the data. Such an analysis gives 

J 0.0107(12) from F| (x 2 = 1-73/2 dof), , , 

° od ~ \ 0.0120(12) from F| { X 2 = 1.18/2 dof). 1 ' 

An overall estimate, obtained by averaging these two values, is 

a od = 0.0114 ±0.0012. (12) 

The data and fits leading to this result are displayed in Figure |2|. 

It is possible to perform an interesting exercise with the data on F\ obtained in || |10 



Recall that this data was obtained through simulations in square lattices where the existence 



of the plateau due to interfaces was not clear. Nevertheless, a rough correction may be 
performed if one argues that the interfaces could be either parallel or perpendicular to the 
z-direction, and hence the net result would be to multiply p™ n in eq. @ by a factor of 
two. Applying this correction, one obtains from the square lattice results a value for a ol i 
surprisingly close to our result quoted in eq. fll^l). Note, however, that in the absence of a 
clear plateau around P™ m this result is accidental. 

In summary, we have performed a measurement of the order- disorder interface tension, 
a d, m the two-dimensional seven-state Potts model and found a value in close agreement 
with the result expected on the basis of perfect wetting. The measurement technique used 
was the so-called histogram method. We have demonstrated that using rectangular lattices 
of shape L x 2L where L ^> £, one may overcome previous objections to this method and 
clearly obtain a plateau at the minimum of the histogram. This, and the use of capillary 
wave corrections to the interface shape are sufficient condition for a reliable application of 
this technique to the measurement of interface tensions. 

We thank Andre Morel for computing the probability distributions for the seven-state 
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of us (SG) would like to acknowledge the hospitality at TFT, Helsinki while this paper was 
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